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Abstract We theoretically study the vortex formation from the collision of the 
domain walls in phase- separated two-component Bose-Einstein condensates. The 
collision process mimics the tachyon condensation for the annihilation of D-brane 
and anti-D-brane in string theory. A pair annihilation leaves the quantized vortices 
with superflow along their core, namely superflowing cosmic strings. It is revealed 
that the line density and the core size of the vortices depend on the initial distance 
between the walls. 
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1 Introduction 

Topological defects or solitons appear not only in condensed matter physics but 
also in cosmology and high-energy physics. It is interesting to connect the physi- 
cal phenomena of topological defects in condensed matter to those in cosmology 
or high-energy physics 1,2 , since the latter are difficult to be realized in experi- 
ments. Atomic-gas Bose-Einstein condensates (BEC) are ideal systems for testing 
the physics of topological defects in theoretical and experimental aspects. One 
of the theoretical advantage is that the topological defects can be accurately de- 
scribed using the mean-field theory. Moreover, the development of the experimen- 
tal technology enables us to control and directly visualize topological defects in 
the atomic clouds. 
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Recently, we proposed that analogues of D-branes 3 , solitons in string theory, 
can be realized in phase- separated two-component BECs 4 , where the two systems 
are connected through the nonlinear sigma model 5,6 . In this model, the roles of D- 
branes and fundamental strings ending on them are played by domain walls and 
quantized vortices ending on the wall, respectively. There exist many advantages 
to study these solitons in this system, since it is possible to investigate experimen- 
tally and theoretically the dynamics of these solitons, such as oscillation modest 
and instability 9,10 of these solitons. For example, Anderson et al. 11 observed cre- 
ation of vortex rings via dynamical (snake) instability of a dark soliton, where a 
nodal plane in one component was filled with the other component forming a pair 
of domain walls. In our context, this experiment can be regarded as a simulation of 
the string nucleation via the D-brane-anti-D-brane annihilation due to the tachyon 
condensation — . 

In this paper, we study the vortex formation by a pair annihilation of domain 
walls of phase-separated two-component BECs in uniform systems. The config- 
uration of the subsequent vortex rings was strongly affected by the geometry of 
the trapping potential in the experiment 11 . Since we are not interested in such a 
finite size effect, we confine ourselves to an uniform system and show that vortex 
formation becomes more complex (see Fig.Q}. The vortex formations are strongly 
affected by the existence of the filling component, which was neglected in the ex- 
periment 11 . The filling component strongly affects on the line density and the core 
size of vortices nucleated after the annihilation. 




Fig. 1 Dynamics of the vortex formations after the collision of the two domain walls. Figures 
show a normalized isosurface of the density difference (n\ — ri2)/n = 0.1. The upper and lower 
figures are the snap shots of the time evolutions from the stationary state with (a-d) = 0.8^1 
and (e-h) \i2 — 0.92/ii, respectively. The darkness of the surface plot refers to the phase S2 of 
the condensate wave function of the component 2, which is trapped along the vortex cores of the 

component 1 [(d) and (h)]. The box size is (51.2 x £) 3 with the healing length ^ = %/ yjmgn\. 
The unit of time is 1/ (Oq = %/ 
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2 D-brane and anti-D-brane in two-component BECs 

The two-component BECs are described by the condensate wave functions (order 
parameters) Wj = y^e^ 7 ' in the mean-field approximation at T = 0, where the 
index j refers to the component j (j = 1,2). The wave functions obey the coupled 
Gross-Pitaevskii (GP) equations^ 



where we use the particle masses mj, the intra-atomic interaction parameters gjj, 
and the inter-atomic interaction parameter gn. When the condition gn > ^gngn 
is satisfied—, the two condensates undergo the phase separation making the in- 
terface layer between domains. In the domain of the component j far from the 
interfaces, we have ttj « = fij/gjj and « (k ^ j). We can regard a plane 
n\ — ri2 in the interface layer as an analogue of D-brane 4 . We use the parameters 
m = m\ = m2 and g = gn = g22 — 0.5gi2 similar to Ref. 4 , where the interface 

thickness is ~ £ = h/ Jmgrfi. 



Let us consider two flat interfaces perpendicular to the z-axis, where a domain 
of the component 2 is between a D-brane at z — —R/2 < and an anti-D-brane at 
z = R/2 > 0. If the distance between the D-brane and the anti-D-brane is large 
enough, the correlation between the two domains of the component 1 in z < —R/2 
and z > R/2 is negligible, and the phase difference dS\ = S+— S- between the 
two domains is arbitrary determined, where we defined the phases as S\ — S- in 
z < -R/2 and Si = S + in z > R/2. 

When the two interfaces are close to each other, the phase difference dS\ = 
S+ — S- becomes important. Since the phase S\ (z) must continuously change in 
—R/2 < z < R/2, the phase difference dS\ ^ yields the gradient energy from 
the spatial variation of S\ (z) after the collision of the two interfaces. If the phases 
—7Z<S\(z)<7Z are mapped to a unit circle, there are two ways from S- to S+, 
clockwise or counterclockwise. The shorter routes are preferred to decrease the 
gradient energy. Especially, for \dS\ — n, there are two possible routes from S- 
to S+. The possibility of different routes can leave topological defects along the 
plane, at which the interfaces collide. This situation is similar to that of the dy- 
namical (snake) instability of nodal planes (dark solitons) for single-component 
BECs in two or three dimensional systems. However, the instability would be 
suppressed for sufficiently large R due to the presence of the component 2, with 
which the nodal plane is filled. 



3 Emergence of excitations with complex frequencies 

To reveal such an effects on the vortex formation, we first investigate the lin- 
ear stability of the dark soliton of the component 1 whose density dip is filled 
with the component 2, which is called dark-bright soliton 15 . Then the condensate 
wave functions in a stationary state may be written as = i// / -(z)exp(— ifijt/h) 
with real functions The excitations are expressed by the perturbations SWj = 




2m 



-v 2 + J Lgj k m 2 vj, 



k 



(1) 
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Wj — from the stationary state. Because of the symmetry of the stationary state 

W9, we have SWj = - tyr* with 9/f(r,t) = uf(z)exp{ikx-imTiVjt/h), 
where fij is the chemical potential of the component j, where we may neglect 
the coordinate y for simplicity in the linear analysis. By linearizing the GP equa- 
tions (Q]) with respect to 5^, we obtain the reduced Bogoliubov-de Gennes (BdG) 
equations 



-gn¥i gnViVi -gn¥2¥i\ 



gn¥t 

\g\2V1V2 -gnWiWi gllWl 



(2) 
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where u u x , w+, w 2 ) r and hf = ^ - +gi^ 2 - lij + ZigjiWi- When 

the frequency has a non-zero imaginary part Im [&>(&)] 7^ 0, the mode is exponen- 
tially amplified. We can obtain W®, SWj and (0 by numerically solving the GP Eqs. 
(Q) and the BdG Eqs. ©. 
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Fig. 2 (a) The density profiles of the phase- separated two-component condensates with two 
interfaces, [(black broken line) n\ and (black solid line) for — 0.80/ii, and (gray broken 
line) n\ and (gray solid line)«2 for ^2 = 0.92^1]. The component 1 has a node and n\ = atz = 0. 

The units of the vertical and horizontal axes are n® = fij/g and ^ = h/ ^mgn®, respectively, (b) 

The dispersion relations of the interface modes in the states of (a). The units of the vertical and 
horizontal axes are coq = \i\/h and 1 , respectively. 



Figures [2] shows the profiles and the dispersion co(k) for typical two val- 
ues of < 1. The particle number of the component 2, thus, the distance R 
between the branes increases with \i2l\i\. For < 1 , the density ri2 of the 
filling component is smaller than the bulk value n^dXz — 0. The imaginary part 
Im [&>(&)] 7^ appears for a range of k and has its maximum at the wave number 
&max- The maximum value monotonically decreases and the dispersion co(k) be- 
comes eventually real for —> 1- In other word, the excitations with complex 
frequencies gradually emerge when the D-brane and the anti-Dbrane approach to 
each other for < 1 . 

The emergence of modes with complex frequencies means that the system is 
dynamically unstable. The unstable modes are exponentially amplified in the lin- 
ear stage of the instability. The amplification of the unstable modes with finite 
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wave numbers k ^ induces superflows, where the atomic interaction energy is 
transduced to the gradient energy. As a result, the instability leads to the vortex for- 
mations as shown Fig.Q] Such a process of the instability due to the unstable mode 
with complex frequency mimics the tachyon condensation caused by the tachyonic 
mode with complex mass m 2 < in string theory 12 . The tachyon condensation can 
leave lower dimensional topological defects after the annihilation of D-brane and 
anti-D-brane. In our case of the phase- separated two-component BECs, the anni- 
hilation of the 2-dimensional defects (domain walls) leaves 1 -dimensional defects 
(quantized vortices). Since the time scale of the amplification of a 'tachyonic' 
mode is given by ~ 1 /Im [(«(&)], the modes with the wave number ~ /c max are 
dominantly amplified, which should be related to the vortex formations. 



4 Vortex formation from the annihilation of the D-brane and the 
anti-D-brane 
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Fig. 3 Schematic diagrams of the vortex formation in the collision process, (a) The planer 
object bounded between the D-brane and the anti-D-brane transforms into a quantized vortices. 
The top figure shows the flat planer object in a stationary state. Black and gray lines show the 
nodal plane with n\ = (inner wall) and D-brane (or anti-D-brane) with n\ = ri2, respectively, 
(b) The reconnections of the rings of the 'half-quantized' vortex living on the planer object. The 
sign of the vortex ring (tunnel) is represented with solid or broken line. The gray background 
refers to the planer object and the white regions represent the tunnels through the object. The 
reconnection between the rings with different signs yields a bridged vortex, (c) Vortex formation 
after the growth of the rings with random distribution. The random distribution causes a distorted 
vortex line formation. 



In the stationary state without the unstable modes, the nodal (n\ = 0) plane 
exists between the branes (n\ = ri2 planes). Then the region with n\ < n>i bounded 
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between the branes forms a planer object [Fig. 13a) top]. When the unstable modes 
are excited and amplified locally, the planer object is broken and there appear 
tunnels through the planer object as shown in Fig. Oa) (middle). The tunnels are 
characterized by the direction of the phase gradient perpendicular to the planner 
object. Since the phase Si changes from to % or — n through the tunnel, the 
'edge' of the tunnel form a pseudo-vortex which mimics a half-quantized vortex. 
Although a pair of the the pseudo- vortices with different circulations annihilates 
and destroys the n\ < ri2 region, a pair of the vortices with the same circulations 
leaves the n\ < ri2 region as a single-quantized vortex whose core is filled with 
the component 2. These processes look like reconnections between the rings of 
the half-quantized vortex, which live on the planer object [Fig. Htb)]. Then the 
reconnections make a larger rings or leave a bridge of single -quantized vortex. 

Since the initial perturbation can be taken to be arbitrary in the uniform sys- 
tems, the positions and the signs of the tunnels would be distributed randomly 
[Fig- Etc) left]. Then, the growing tunnels form the mesh structures [Fig. [3](c) 
(center)], which causes a distorted vortex line formation [Fig. 0(c) (right)]. 

Now, we can estimate the vortex line density n v after the collision. Since the 
unstable modes with k ~ kjYisLx <ire dominantly amplified, the mean distance be- 
tween the neighboring vortex ring would be of order ~ 1 /k mSLX in the early stage 
of the instability. As a result, we have the vortex line density n v ~ & max per unit 
area along the z — plane. Since & max is decreased with the resultant vortex 

line density increases with 

These interpretations explain well the numerical results in Fig. Q] The numer- 
ical simulation was done by solving the GP Eqs. (Q]) from the stationary states 
for = 0.8, 0.92 with a random seed for the initial perturbation in a periodic 

system along the x-y plane. Since the particle number is conserved in these simula- 
tions, the particle density of the component 2 per unit length along the vortex line 
becomes larger for weak instability with small k m£iX . The fluctuation of the phase 
S2 of the component 2 causes the superflow along the vortex line. Such a vortex 
may be called as a core-flow vortex or a superflowing cosmic string, analogue 
of the superconducting cosmic strings 16 . For = 0-92 in Figs. Q] (bottom), 

the particle number of the component 2 is so large that we cannot define well the 
vortex core in contrast to the case of — 0.80 in Figs. Q] (top). 

In the experiment by Anderson et al.—, the snake instability was induced by 
removing the filling component selectively. This experimental technique may con- 
trol the relative velocity between the D-brane and the anti-D-brane. If they ap- 
proach to each other slowly, the dynamical instability emerges gradually and thus 
the vortex formations are suppressed. The vortex line density n v after the annihi- 
lation increases with the relative velocity. Roughly speaking, the density n v takes 
the maximum value ~ fc max (jU2 = 0) ~ when the relative velocity is larger 
than ~ £ /Tmin, where T^n = l/Im[o)(& max )] is the growth time of the instability 
for = and t, is the thickness of the interface layers. 



5 Conclusion 

We have studied the vortex formation after the pair annihilation of the domain 
wall and the anti-domain wall in phase- separated two-component BECs in uni- 
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form systems. The vortex line density after the annihilation is qualitatively es- 
timated, which may reach to its maximum value for a sufficiently large relative 
velocity between the walls. The nucleated vortices after the annihilations have 
superflows along their core regions. These phenomena can be experimentally re- 
alized in phase-separated condensates whose size along the walls is larger than 
the length scale ~ l/k mSLX . The detailed analysis of this study will be reported 
elsewhere soon. 

It is interesting to mention about the relation to the defect creations by collid- 
ing interfaces in the experiments of superfluid 3 He-ii, where A-phase domain is 
between two B-phase domains forming the A-B and B-A interfaces as analogues 
of branes. The underlying mechanisms of vortex formations proposed in this work 
might be applied to the vortex formation in the 3 He experiments. In our case of 
two-component BEC system, vortex formations are essentially triggered by mak- 
ing tunnels in an unstable soliton of the 1st component hidden in the domain of 
the 2nd component. Since such an unstable soliton is possible in 3 He-B 2 , we may 
expect that the similar vortex formations occur in the pair annihilations of the A-B 
and B-A interfaces in the 3 He experiments if the soliton is hidden in the A-phase 
domain. The theoretical and numerical analysis in this direction are the interesting 
future works. 
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